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I. 


INTRODUCTION 


Solution  of  complicated  boundary  value  problems  in  the 
spectral  domain  is  a  very  well  known  technique.  Several  years  ago 
Itoh  and  Mittra  used  spectral  analysis  to  solve  planar  geometry 
problems  occurring  in  microwave  integrated  circuits,  [1] .  There 
have  been  numerous  papers  subsequently  published  utilizing  this 
method,  [2].  In  this  paper,  Fourier  analysis  is  used  to  formulate 
the  solution  for  reflected  and  transmitted  surface  current 
densities  in  back  to  back  shielded  microstrip  structures  coupled 
by  a  transverse  slot  in  their  common  wall.  By  applying  Galerkin's 
procedure  in  the  spectral  domain  on  complementary  field 
quantities,  a  solvable  system  of  equations  containing  the  desired 
reflection  and  transmission  coefficient  results. 

II.  THEORY 

Application  of  the  Galerkin  method  on  the  spectral  domain 
electromagnetic  quantities  supported  by  the  structure  illustrated 
in  Figure  1  requires  a  matrix  equation  which  expresses  electric 
and  magnetic  fields  in  terms  of  complement  source  densities.  The 
appropriate  expression  is  given  in  equation  (1) . 


(la)  0 

E2  <a,ha,P) 

Gu(ct,P)  G12(o,p)  G13  (a,  P) 

k  (a, h  ,  P) 
z»  a 

~  (la)  _ 

Hz  (a,  o,p) 

3K 

G21<a,P)  G22<a,p)  g23«x,P) 

(a,o,p) 

~  (lb)  a 

Ez  <a,-hb,p) 

G31(a,  p)  c32(a,P)  g33«x,P) 

(1) 


1 


Explicit  expressions  for  the  G^fa,  (J)  ,  along  with  an  outline 

of  their  derivation  are  given  in  Appendix  A.  In  this  derivation, 
we  have  assumed  narrow  microstrip  lines  and  a  narrow  coupling 
slot.  These  acceptable  restrictions  allow  us  to  assume 
uni-directional  surface  current  densities  which  have  a  Maxwellian 
distribution  in  the  direction  of  their  narrow  dimension,  [3] . 
Recent  theoretical  and  experimental  research  by  Dunleavy  has  shown 
that  the  narrow  microstrip  line  current  variation  mentioned  above 
is  very  accurate,  [4] .  Unfortunately,  no  such  verification  exists 
for  a  magnetic  current  on  a  narrow  aperture.  However,  since  a 
narrow  microstrip  line  supporting  an  electric  current  and  a  narrow 
aperture  supporting  a  magnetic  current  are  eseentially  dual 
quantities,  [5],  it  seems  that  an  assumed  Maxwellian  distribution 
for  the  aperture  problem  would  be  accurate. 


III.  SPECTRAL  DOMAIN  REPRESENTATION  OF  CURRENT  DENSITIES 

The  detailed  geometry  of  this  problem  is  given  in  Figure 
2.  In  the  spatial  domain  the  two  electric  surface  current 
densities  away  from  the  discontinuities  are  written  as: 


K  (x',h  ,*') 

Z A  * 

VX,'"VZ,) 


-  Te  gb(x’) 


ga(x-) 


z  *>z 

O 


z  '<>-z 

o 


where 


(2) 


2 


g^x')  = 


KW, 


{i-  [^Vf'2 


w. 


Otherwise 


wi  wi 

—  <  x’-x  <  — 

2  1  2 


(i  =  a,b) 


(3) 


It  should  be  mentioned  that  y  and  y  are  known  constants 

za  zb 

for  their  respective  guiding  structures.  Their  values  may  be 
obtained  from  the  literature,  [2] .  In  the  vicinity  of  the 
discontinuities,  all  surface  current  densities  in  equation  (1) 
must  be  expressed  in  terms  of  appropriate  basis  functions. 


K  (x’,h  ,z')  and  K  (x',-h.,z')  are  written  in  the  spatial 

za  a  zb  d 


domain  as 


K 


N 

(x* ,hi, z' )  =  X  \  9i 

n*l 


(x1 )  f  ( z  * )  i  =  a,b 


(4) 


Where  the  I  are  unknown  amplitudes  and 


f  <z') 

nt 


sin  [  K(z ' ~zn-l)  ] 
sin  (Klz) 

sin  [  K(zntl~z 1 )  ] 
sin(Klz) 

0 


n-l 


<  Z'  £  Z 


Z  £  Z  '  <i  z 

n 

Otherwise 


n+1 


(5a) 


For  n  *  1 
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(  sin  [  K  (±zrl-z')] 


•  i  =  \ 


f.  <z')  = 


o  z 


sin  (Kl?) 


(±«  1  "  lb)  ' 


±z  T1  £z  '£±z 

o  z  o 


Otherwise 


(5b) 


where  1  =  Z  -Z  -  Z  -Z  ,  .  The  magnetic  current  on  the  coupling 

z  n+1  n  n  n-l  3 

aperture  is  represented  as 


Mx  (x' / 0/  z ' ) 


i 

n-l 


v  f  (x')  g.(z') 

n  n  o 


(6) 


where  fn(x')  and  g0(z')  are  given  as  equations  (7)  and  (8), 
respectively,  and  the  Vm  are  unknown  amplitude  coefficients 


sin  [  K(x’-xn_x)  ] 


sin(Kl  ) 


fn (X '  ) 


=  |  sin  [  K(xn+1-x')  ] 


sin(Kl  ) 

X 


X  ,  <1  x'  £  X 
n-l  n 


x  S  x'  S  x  , 

n  n+1 


Otherwise 


(7) 


(  \  -1/2 


-  -  £  z'  £  — 
2  2 


(8) 


Otherwise 


4 


We  also  see  that  lx  =  xn+1  -  xn  =  xn  -  xn_!  and 

xn  =  (n-1)  lx.  Under  these  conditions  the  requirement  that  the 

magnetic  current  is  non-zero  only  on  the  coupling  aperture  is 
satisfied.  The  next  step  in  the  analysis  is  to  obtain  the  Fourier 
transforms  of  (4)  and  (6) .  The  transforms  of  interest  are  written 
as : 


V°'h>  ■  JJ  vx,'>v 


*•)  e)OX’  e1""  dx’dz' 


(9) 


KiB(o,-h;-P>  -  JJ 


jax'  jBz* 

z')  e  e  dx'dz' 


(10) 


Mxja,0,ft)  =  JJ  <x',0,z') 


jax'  jPz' 

e  e  dx'dz' 


(11) 


For  K  (a,h  ,P)  we  obtain 

z_  a 


V“-VP>  -  J  g,(x 


jax1 

' )  e  dx ' 


[  t I 

n“2 


jBz  • 

f  (z')  e  dz' 

n 


z' 


J  Ix  fx(z')  ejPz'  dz'  +  J  (e1*  -  Re  K  )  e'r~  dz 


zo~lt 


V»>vp)  - 1*  [i;  + 1;  + 1;] 


(12) 


2 

7CW 


.  [2^1  C 


jo*'  ,  , 
e  dx* 


(13) 
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m  zr»  zn+l  , 

^  f  sin  K(z'-2n-l)  jPz,  r  sin  K(zn+1-z')  jPz. 

=  7,1  I  - ; -  e  dz'  +  I  - ; - - -  e  dz 

^  na  J  sin(Kl  )  J  sm(Kl  ) 

n-2  ,  z  ,  z 

zrr  1  zn 


■;  ■  I  ' 


jp)z<  <-r“ 
-  Re 


+  jP)Z' 


After  some  algebra  I  ,  I2/  and  13  are  found  to  be: 


a  3°*.  /  wa 

X1  =  e  J4a  T 


(14a) 


u 

-  /  ^  £> 


jBz' 

(z  ' )  e  K  dz  • 


(14b) 


a  j2&  V  J 

12  =  ^  ^  6 

(K  -p  )  n-2 


(17a) 


JP<z0_1x> 

zl  e 

.a  _ 

(K2-p2)  sin(Klz) 


[  IK 

K-e 


(j^sinKl  -  kcosKl  ) 


(17b) 


To  evaluate  IA,  some  special  considerations  are  necessary. 


I4  is  written  alternatively  as 
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-  J [<' 


iP>  z '  <-y" 

-  Re  ' 


*  +  jP)  Z' 

a 


and  we  assume  a  lossless  system  so  that  YmS  is  a  purely  imaginary 

za 

,  ms  ms  , 

quantity,  y  =  j  k  .  The  Fourier  transform  of  the  product  of  two 

za  za 

functions  is  equal  to  the  convolution  of  the  individual  Fourier 
transforms  of  the  two  functions,  [6] .  Using  this  idea  to  evaluate 
(18)  shows 


ms 

ms  \ 

\  J 

r  z' 

-J*,  T  V 

}  * 

Le  -  Re  *J  J 

op  •• 

J  O(z'-Zo)  e^dz1  “  ^(P)  +  “ 
j  e**1’ e^dz-  -  2*8  (P  ♦  \) 


1  iPz0 


00  ms  /  \ 

_R  J  e  ^  e^’dz’  =  -27CR5  (P  -  <j 
— «• 

,  a 

From  equations  (19)  and  (20)  we  express  I4  as 


(20a) 


(20b) 


(20c) 
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<  - 2*  J  e)lI“  I***’ +  * 


S  (fi-X  +  k“)  -  r5  (p-X  -  k™) 


dX 


=  2n 


r  i  jaz0 

f  L1®5  »<»•>  8  (p-*  +  *"')  +  - - S^-Ji  +  k”) 

J  za  z* 


jX 

j*-zQ 


1  J  0  1 

-TCRe'  Z°  8(X)  6  (p-X  -  k"s)  -  — -  8  (p-X  -  k”s)  J 


jX 


dX 


I!  -  2*2  s(p  +  k“)  -  2k2  r8  (p  -  o 

i(K')  . 


+  2n 


-e 


j  (p+k“)  *0 


+  R 


j(p<)  j  (p-C) 


(21) 


Substituting  (16),  (17),  and  (21)  into  (12)  shows: 


w 


jOxa  -j(p+fc™)z0 

K  <“'h,-P>  -  -  *  6  .  e - — - J  |  a  ^-1 

i(K')  1  2 


i«*.  -j  <M“)  , 

2n  e  e 

+  R  - - - - -  J 


w\  *  Jo  (°  -f)  J, 

—  L 


j  (p-C) 


3PZn 


a  —I  + 

2 


(k2-P2) 


n=2 


I  e 

n. 


jox,  jp<z0-it) 
e  e 

+  i.  - - -  j 


la  (k2-P2)  sin(klz) 

2jc2  S  (p+k”)  -  r5  ( P~k™* )  j  ej<XX‘  Jo 


Wa  jpiz 

a  “  K  -  e  1  (jPsinKl  -  kcosK  1  ) 


w 


(22) 
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In  a  completely  analogous  manner  K  (a, -h.,B)  is  found  to  be 

zb  b 


k  (a, -h  ,  P)  =  -t 


2k  e 


jaxb  j  (  P+k”  )  zg  (  w  ^ 


j  (p<) 


a 


-:iP<zo-V 

I.  e  e 

J 

a  — 

0 

l  2j 

(k2-P2)  sin(kl  ) 


K  +  e  (j^sinKl  +  kcosKl  ) 


jaxh 


jB  e  J  a 


2  ^aXb 

I  e  +  2i  T  e  J 

2  n2\  c 


(k  -p  )  - 


w. 


a  y]  6  (  |3<s  ) 


(23) 

The  discretized  microstrip  currents  near  the  discontinuity 

region  are  divided  in  'N'  partitions  in  both  guides  'a'  and  ' b' . 

More  severe  disruption  of  the  current  occurs  in  guide  'a'  due  to 
the  open  end  upon  which  Kz  is  incident.  If  'N'  partitions  are 

required  to  accurately  characterize  Kz  then  the  same  number  will 

represent  K  in  its  respective  discontinuity  region. 

zb 

The  Fourier  transform  of  M  (x',0,z')  is  expressed  from  (6), 

xa 

(7) ,  (8)  and  (11)  as 


m  (a,  o,p) 


fn(x’)  g0(z’) 


jax'  jpz1 

e  e 


dx '  dz ' 


(24) 
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and  clearly 


t/2  /  .i2V1/2 

=  J  —  ( -  >  [cos^z*  +  jsin($z']  dz' 

-f/7  V  U  J  / 


-J=r 


cospz'  dz' 


2z '  -t 

^  =  cos<|>  =>  dz '  =  —  sin<(>  d<|) 


4  =  f  — - —  cosp  ~  cos<|>  -  sin<j> 

J  JCt  sin<j>  v z  J  2 


!<  -  ^  J  cosP(|-  cos*]  d*  =  Jo[P  ^ 


M 

I5  =  X  V  f  f 
5  n  J  n 


jOx' 

(x ' )  e  dx ' 


Substituting  for  fn(x')  from  equation  (7)  and  evaluating  the 
result  yields 


m 

i  .  y  v  ^2a  e 

5  ^  n  22 

n-1  (K  -a  ) 
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Finally  then,  M  (a,  0,|3)  is  written  as: 

X_ 


M  = 

xa 

Equations 


j2a  Jo(p  |) 

2  ? 

(K  -a2) 

(22),  (23), 


n-l 


jOXn 

e 


(28) 


and  (28)  are  very  important  quantities. 


IV.  FORMULATION  OF  THE  MATRIX  EQUATION 

The  results  obtained  in  the  previous  section  will  allow  us 
to  solve  equation  (1) .  In  the  following  discussion,  spatial 
and/or  spectral  dependence  of  functions  will  be  implied  and  not 
explicitly  written.  If  there  is  any  ambiguity,  it  will  be  clearly 
explained  from  equation  (1)  : 


~  (la) 

E 

Z 


G  K  +  G10  M 

11  za  12  xa 


H<la’ 

Z 


G,_  K  +  G,_  M 

21  za  22  x. 


~  (lb) 

E 

z 


G„  M  +  G„  K 

32  x,  33  z„ 


Since,  from  Appendix  A  we  know  that  G13  =  G23  =  G31 


(29) 

(30) 

(31) 


0. 


Also,  note  that  the  surface  current  densities  are  expressed  in 

terms  of  summations  of  appropriate  basis  functions  as  derived 

earlier.  Before  going  any  further  with  the  analysis  it  will  be 

convenient  to  write  K  ,  K  ,  and  M  as: 

V  Zb  X3 


(32) 


N 

~a  ~a  ~n 
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K  =  -  T  k  +  /-I  I  K 

zb  1  ~  nt,  zb 

n»l 


M  *  "V  V  Mn 
*,  jL^  n  xa 


where  kQ,  k^,  and  Kz  are  easily  deduced  from  (22) 


~b  ~  n  ~  r 

k  ,  K  and  M 

1'  X 


are  similarly  obtained,  respectively,  from  (23)  and  (28) .  We  next 
multiply  equation  (29)  by  Kq  and  integrate  from  -«>  to  <»  with 


respect  to  a  and  P  for  different  values  of  q.  This  shows: 


JJ  S’  i;u’do  dp  -  JJ  S’  Gu  K  da  dP 


JJ  Kz  G12  M^da  df  q  =  1,2,  ...,N 


JJ  Kq  Ezla)da  dp  =  0 


Because,  from  Parseval's  theorem 
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12 


~q  ~  <la) 

Since  it  is  clear  that  Kz  and  Ez  are  complementary  in  the 
x-direction.  Using  (32)  and  (34),  equation  (35)  becomes: 
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(36) 


q  =  1,2, . . .,N 


Performing  similar  operations  on  equations  (30)  and  (31)  yields: 
From  (30) 
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From  (31) 
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(38) 

N  00 
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ri”l 

Equations  (36),  (37),  and  (38)  represent  an  inhomogeneous 

system  of  2N+M  equations  containing  2N+M+2  unknowns.  We  get  two 
additional  relations  by  the  requirement  of  electric  surface 

current  continuity  at  (x=ha,  z=*z0)  and  (x=-hb,  z--zQ)  .  Look  at 
Figure  2 . 

From  equations  (2),  (3),  (4),  and  (5)  we  see 


I 


(39a) 
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T  e 


(39b) 


Equations  (39)  represent  the  necessary  relations  for 
solving  for  R  and  T.  We  shall  now  express  (36),  (37),  and  (38) 

more  compactly: 

(36)  becomes 
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$  -  R  Q  +  7-1  L  +  7,  V  F  q  =  1, 2,  .  .  .  ,M 

”  "  n»  ”n 

n-1  n-1 


where 


-R 


g21  k,  dadP 


-  JJ  G„  i‘  dodP 


G„  K,  dad? 


JJ  «x,  5„  \  d0dP 


(38)  becomes 
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-  jj  \  G33  *1  d“dP 

—99 

99 

vq»  *  JJ  G»  \  dod&  <42» 

—OP 

oo 

X,„  -  JJ  Kb  532  \  dddf> 


V.  DISCUSSION  or  NUMERICAL  SOLUTION 

It  is  clear  from  the  results  of  the  last  section  that  the 
problem  under  consideration  is  very  difficult.  Since  no  numerical 
results  will  be  obtained,  a  discussion  of  how  to  do  this  is  in 
order.  Equations  (40),  (41),  and  (42)  contain  integrals  which 

must  be  evaluated  from  to  with  respect  to  two  variables:  a 

and  p.  A  brief  glance  at  Appendix  A  indicates  that  numerical 

integration  of  these  terms  will  surely  be  required.  Although  it 

is  not  obvious,  these  integrals  vary  at  worst  as  — — —  ,  —  .  — - —  .  — 

a  a  P  P 

in  one  part  or  another.  The  squared  terms  will  converge  rapidly 
and  the  same  can  be  said  for  the  — ~  and  terms.  Ultimately,  one 

3/2  a  3/2  J 

a  P 

would  need  to  program  a  computer  to  do  these  integrations  to  find 
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how  the  convergence  progresses.  However,  in  [7],  Mittra  and  Itoh 
discuss  so  called  'numerically  efficient  techniques'  for  solving 
boundary  value  problems.  The  approach  used  in  this  report  closely 
resembles  the  algorithms  discussed  there.  As  a  result,  it  would 
seem  that  satisfactory  convergence  is  obtainable  for  the  equations 
derived  here.  Besides  having  to  evaluate  the  integrals  mentioned, 
we  must  also  solve  a  linear  system  of  equations.  There  are  many 
viable  alternatives  for  this  part  of  the  problem  and  no  further 
discussion  is  necessary. 

VI .  CONCLUSION 

The  solution  for  reflection  and  transmission  coefficients 
in  slot-coupled  microstrip  lines  has  been  formulated  in  the 
spectral  domain.  Although  no  numerical  results  have  been 
obtained,  pretinent  numerical  considerations  have  been  discussed. 
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Appendix  A 


Derivation  of  the  Green’s  matrix 
Hybrid  mode  analysis  with 

E  (x,  y,  z)  =  -j£OA  +  — - —  VV  •  A  +  —  V  x  F 

jo)ji£  e 


(A-l) 


H(x,y,z)  =  jCOF  - 


jCOjie 


VV  •  F  +  —  V  x  A 


Choose  A  =  zA  and  F  = 


zF  where 

Z 


(x,y,  z) 


(A-2) 


2  IT  (a,y,P)  e  jax  e  jPz  dadfi 

(270  " 


~(i) 

A 

z 


(DUE  |  ^  (i) 

— ]  ¥a  (a,y,p)  (i  =  la, lb, 2a, 2b) 


(i) 

F  (x,  y,  z)  = 


1  ff  ~(1)  n  -jax  -jPz  _ 

- 2  JJ  F*  (a,y,P)  e  e  dadP 


(270  4.' 


(A-3) 


~  (i) 
F 

z 


(^)  (a'y'P>  <i  =  la,  lb,  2a,  2b) 


Aza)(x,y,  z)  and  F^l)(x,y,z)  must  satisfy 


(A-4) 


(v’ +  *’) 


A*1’  (x,  y,  z) 


(i) 

F  (x,y, z) 
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Substituting  (A-2)  and  (A-3)  into  (A-4)  implies: 


Appropriate  solutions  to  (A-5)  in  each  region  of  the  structure 
shown  in  Figure  1  are: 

=  AUa)  sinh(yiay)  +  B<la>  cosh(auy) 

-  5“’  sinh[irjr*>)] 

(A-6)  ijrj[lb)  =  A<lb>  sinh(yiby)  +  B<lb>  cosh(aifey) 

=  i'2"’  sinh^M] 

Vpla)  =  C<la>  sinh(yiay)  +  D(la>  cosh(aiay) 

Vp2a)  =  N(2a>  cosh[y2a(y-b)] 

(A-7)  =  C<lb)  sinh(Ylby)  +  D<lb>  cosh(aiby) 

-  N(2b)  cosh[y2b(y+b)] 

Note:  The  ~  notation  means  Fourier  transform  domain. 
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( 


The  transformed  fields  must  satisfy 


(A-8 ) 


i(lq) 

z 

(a, 

=  El2q) 

z 

(a,hq, 

P> 

~  (lq) 
E 

X 

(a, 

~(2q) 
=  E 

X 

(a,hq> 

P> 

~  (lq) 
H 

z 

(a, 

=  H<2q> 

z 

(a,hq, 

P) 

-aq) 

X 

(a, 

VP> 

~(2q) 

-  H 

X 

«x,hq 

,p)  =  KZq  (a,hq,p) 

E(lq) 

X 

(a, 

O.P) 

=  0 

i(lq) 

Z 

(a, 

0,P) 

=  M,<a, 

0,p> 

(q  =  a,b) 

K  (a,h  ,B)  is  the  Fourier  transform  of  the  electric  surface  current 
zq  q'K 

on  the  qth  microstrip.  Mx«X,0,p)  is  the  Fourier  transform  of 

magnetic  surface  current  on  the  aperture  which  is  common  to  the 
two  microstrip  structures. 

At  this  point  in  the  derivation  it  should  be  noted  that 
symmetry  about  the  line  x=0,  y=0  exists  between  guides  (a)  and 

(b)  .  In  order  to  derive  Gn(a,j$)  and  G12(a,|J),  we  apply  the 
boundary  conditions  (A-8)  and  express  E<U)  (a,h  ,  ($)  appropriately. 

Z  a 

G32(a*P)  and  G33(a,P)  are  obtained  from,  respectively,  G12(a,P) 
and  Gu  (a,  f$)  as 

G„(a,fJ)  *  G13(a,  p)  |  ,  _h^;  ^ 


(A- 9) 


G33(a,P)  =  Gn(a,P)  |  h> 


-hb;  Kla  -  ku 
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Consequently,  solving  for  G  2  (a,  P)  and  Gu(a,P)  yields  G32  (a,  P) 

and  G33 (a, p)  in  a  direct  manner.  Imposing  the  boundary  conditions 

(A-8)  and  appropriately  manipulating  terms  yields  the  components 
of  the  Green's  matrix  as: 

GU«X,|3)  =  -  jCOM-  (*L-  P2)  (k2a-P2)/  {a2p2(kla‘  k2a) 

-  ^al^a-  32)  "  rLkia(k 2a" 

+  YiaY2a(kL"  P2)  (k2a“  P2)  {  kL  tanh[y2a(h-b)]  coth(yiaha) 

+  k2a  tanh(yiaha)  coth[y2a(ha_b)]  }  } 


g21«x,P)  = 


-jap(k2a"  kla)  COSh(yiaha)  Gn<(X,p) 


°^Y2a(kxa-  P2)  tanh[y2a(Vb)]  -  yia(k2a~  P2)  tanh(yiaha) 
- 


G  =0 

31 


GJ2(a,p)  -  cosh(rla\)  -  2p(k',-  p2)  (k2a-  p2)  (k2a-  k2a)  -f 


fa(a,p,ha) 


(a,p,ha) 


f  (a,P,ha)  _  sin(rl4»a)  {  tanh[r2a(hrb)]  !p2Y2a(ki.~  (kL+  yL) 

-  2,(kL-  kL)  ♦  «2y  Ms  t  L)  (kL-  e2)  -  kLv  Ms  ] 

+  tanh(rlaha)  [k[.r2.r Ms  p!)  '  T'-)  P2) 
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-  “2y  „(kL+  1  L)  R-  P2)  ]  -  k2aYi,Y2,(kL  -  P2)  =°th[v2ih.-b)]  } 
+  coshRhJ  [ct^R-  kjj  -  Y1J2,k2.(ki.'  P2) 

5.(°<P'h.)  *  “P2(kl,-  k2.)  +  k22.Y2,(kL'  P2)  +  kLY2,(k2.-  P 

+  YuY;,k!a(kL-  P2)  R-  P2)  tanh(Yla»J  coRjV13)] 

-  Y  1aY2akia(k2a-  P2)  R*  P2)  coth(Y1,hJ)  tanRjV13) 


G22(a,p)  =  fj(a 


,p)  /  t(a,p> 


where 


Tl(a,p) 


sinh 


Rh.) 


oa\iy 


)  -  R+  y2,)(“2+  p2) 


p2  kL- k2.)  -  fc2.+ y2 


,  2  2 
+  k  Y. 

1  •  0  a 


la 

we,  Yr-  ^ 


la#^""  la 


R-  p2) 

COSh(Yiaha)  COth  Rh.)  }  *  {  ap(kja  -  k2J  / 


R- 

R2.R  -  P2)  “th[Y2,(h._b) 


-  ei .*,«■ 


R- p2)  coth(Ti»h*)]  } 


+  jog 


R-  p2) 


,fT  ,)  T*  ®lRh.) 

0  1  lAJ "  U  7\ 

Yl,(.k2,'  P  j 


tanh 
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t  O  «  UJ 


-  T  • 

R-  »2) 

(m  [ 
1  H 

y2{klr  P2)  tanh[rjh.  b) 

“  Yla(kL"  P2)  tanh(Ylaha)  ] 

-  &  R-  kL)  /  U,t2,R 

coch[v2,(h*  b)] 

-  E  V 

laha 

R  ■  P2)  cothRh,)  ]  } 

32(a,P)  is  obtained  through  equation  A-9. 

G13  (a,  (J)  »  0 
G23  (a,  f})  =  0 

(a,^)  is  obtained  through  equation  A-9. 


24 


